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Linear stability and the non-modal transient energy growth in compressible plane Couette 



^ ■ flow are investigated for two prototype mean flows: (a) the uniform shear flow with con- 

. stant viscosity, and (b) the non-uniform shear flow with stratified viscosity. Both mean 

O 

I flows are linearly unstable for a range of supersonic Mach numbers (M). For a given 

(— I ' M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow 

than its non-uniform shear counterpart; for a given Re, the dominant instability (over all 

, stream-wise wavenumbers, a) of each mean flow belongs different modes for a range of 

' supersonic M. An analysis of perturbation energy reveals that the instability is primarily 
O ' 

Q ■ caused by an excess transfer of energy from mean-flow to perturbations. It is shown that the 



energy-transfer from mean-flow occurs close to the moving top-wall for "mode I" instability, 
whereas it occurs in the bulk of the flow domain for "mode 11". For the non-modal transient 
growth analysis, it is shown that the maximum temporal amplification of perturbation en- 
ergy, Gmax> and the corresponding time-scale are significantly larger for the uniform shear 
case compared to those for its non-uniform counterpart. For a = 0, the linear stability op- 
erator can be partitioned into £ ~ £ + Re^Cp, and the i?e-dependent operator Cp is shown 
to have a negligibly small contribution to perturbation energy which is responsible for the 
validity of the well-known quadratic-scaling law in uniform shear flow: G{t/Re) ~ Re^. 
In contrast, the dominance of Cp is responsible for the invaUdity of this scaling-law in non- 
uniform shear flow. A reduced inviscid model, based on EUingsen-Palm-type solution, has 
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been shown to capture all salient features of transient energy growth of full viscous problem. 
For both modal and non-modal instability, it is shown that the viscosity-stratification of the 
underlying mean flow would lead to a delayed transition in compressible Couette flow. 

PACS numbers: 47.15.Fe, 47.20.Ft, 47.40.Ki 



I. INTRODUCTION 



The transition to turbulence has its genesis to the growth of small disturbances in an other- 
wise undisturbed laminar flow. Hence, an understanding of different mechanisms of instability 
growth is important to determine the transition scenarios that eventually lead to turbulence in flu- 
ids. The linear stability theory, via the standard normal-mode approach, is the starting point of 
such analyzes which predict the onset of the growth of a small perturbation imposed on a laminar 
flow [jl]. However, there are flow configurations that are stable according to the linear stability 
theory (i.e., the critical Reynolds number is infinity, Rccr = oo) for which the experiments show 
a finite i?ecr(<< oo). It is reasonable to assume that such sub-critical flows may be destabilized 
by the nonlinear effects that are neglected in the linear theory. Is there any linear mechanism 
that causes an infinitesimally small perturbation already present in the flow to grow substantially 
for a short time? If this is true, then the nonlinearities could take over subsequently to trigger a 
flow-transition. 

Indeed, following the seminal work of Orr [|20, it has recently been realized that one should 
investigate the short-time dynamics of the flow without a priori assuming the exponential time- 
dependence for perturbations. The key idea is that even though each eigenmode decays in the 
asymptotic limit (t oo) for a stable flow, a superposition of such stable eigenmodes has poten- 
tial for large transient enemy growth before t hey can be stabilized by the viscosity. Such transient 



growth analyzes y, 



alii 



1211 have revealed that a flow can sustain large am- 



plification of perturbation energy even if the flow is linearly stable. In mathematical terms, the 
underlying linear stability operator is non-normal (non-self-adjoint) [|4|, |5i [l3|] which is responsi- 
ble for transient energy growth. This is a possible route to flow transition for subcritical flows 



IJ. 



which has become an active field of research during the last ten years 

It is known that small changes in the mean-flow can be stabilizing or destabilizing which is 
an attractive avenue from the viewpoint of controlling or manipulating instabilities. A recent 



work ll^ clearly underscores the effects of mean-flow variation on the stability of incompress- 
ible plane Couette flow- using the concept of pseudo-spectra Isl, 0], these authors showed that 
relatively small changes in the mean flow could be destabilizing. Another important issue in 
stability research is the possible role of viscosity stratification on instabilities which has a stabi- 
lizing effect, leading to a delayed transition. In incompressible non-Newtonian fluids, the role of 
viscosity-stratification in delaying transition is currently debated for which we refer to a recent 



work nnn . 

The above issues have not been investigated for compressible fluids in a systematic manner. 
In this paper, the linear stability characteristics and the non-modal transient energy growth in the 
compressible plane Couette flow are analyzed for two prototype model problems: (a) uniform 
shear flow with constant viscosity, and (b) nonuniform shear flow with stratified viscosity. The 
first goal of the present work is to understand the similarities and differences of the modal and 
nonmodal stability characteristics between these two closely related mean flows of a compressible 
fluid. The second goal is to reveal the role of viscosity-stratification on instabilities in a com- 
pressible fluid since we have two prototype mean-flow configurations in which one has a stratified 
viscosity across the channel and the other has a constant viscosity. The third goal is to characterize 
the underlying instability mechanism in compressible Couette flow via an energy analysis. 

This paper is organized as follows. The governing equations and the mean flow are briefly 
described in Section II. The linear stability problem is formulated in Section III, and the related 
results are presented in Section III. The results on the transient energy growth are presented in 
Section IV. The summary and conclusions are provided in Section V. 



II. EQUATIONS OF MOTION AND MEAN FLOW 



Consider a perfect gas of density p* and temperature T* between two walls that are separated by 
a distance h*: the top wall moves with a velocity Ui and the lower wall is stationary, with the top- 
wall temperature being maintained at T/; here the superscript * denotes dimensional quantities, 
and the subscript 1 refer to the quantities at the top wall. Let u*, v* and w* be the velocity 
components in the streamwise (x*), wall-normal (y*) and spanwise (z*) directions, respectively. 
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The conservation equations for the mass, momentum and energy, in dimensionless form, are: 

I = -V.(pu), (1) 



Dui 1 dp 1 



d 

pSj\i + ^(AV • u) 
d c}w. 

+ ^d^-^^ ' + ^^^^ ' ^^""'^ + ^^^^ ' d^- 



(2) 

P^=(1-7)pV.u+^V.(^Vt)+$ (3) 



with D/Dt — {d/dt + u • V) being the material derivative, $ the dissipative shear work, and the 
equation of state is that of a perfect gas: p — pT. We have used the separation between the two 
walls h* as the length scale, the top wall velocity, Ui, and temperature, T^, as the velocity and 
temperature scale, respectively, and the inverse of the overall shear rate, U^/ h*, as the time scale. 
The nondimensional control parameters are the Reynolds number Re, the Prandtl number a and 
the Mach number M: 

p*mh* pi*ci m 

Re = ^ , a = ^ and M = , ^ . (4) 

Here p* is the shear viscosity, C, the bulk viscosity, k* the thermal conductivity, 7 = c*/c* the ratio 
of specific heats, R the universal gas constant and X — ( — 2ij,/3. The bulk viscosity is assumed 
to be zero (i.e., C = 0) such that A = —2p/3 (Stokes' assumption). For all calculations below, 
a = 0.72 and7 = 1.4. 



A. Constant viscosity: Uniform shear flow 

For the unidirectional steady and fully developed mean flow, the continuity and the z- 
momentum equations are trivially satisfied. From the y-momentum equation, it is straightforward 
to verify that the pressure, po — po{y)To{y), is a constant, which is normalized such that po — 1. 
(The subscript is used to designate the mean flow quantities.) The boundary conditions on the 

stream- wise velocity Uo{y) and temperature To(y) are 

C/o(0) = 0, C/o(l) = l, To{0)^T^, To(l) = l, (5) 

with being the temperature of the lower wall. 

For the constant viscosity model (po — constant), the stream-wise velocity varies linearly with 

y- 

Uo{y) = y, (6) 
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i.e., the shear-rate is uniform. It is straightforward to verify that the temperature varies quadrati- 
cally with y: 



To{y) = Tr 

where is the recovery temperature, 



r + (1 — r)y — I 1 



y 



Tr = l + 



(7 - 1)(tM2 



(7) 



(8) 



and r = T^u/Tr the temperature ratio. Note that r = 1 corresponds to an adiabatic lower wall. 



B. Viscosity stratification: Nonuniform shear flow 

For a temperature-dependent viscosity model, for example, with Sutherland's law, 

the streamwise velocity has a non-uniform shear rate. For this case, the mean flow equations, 

with boundary conditions © have been solved numerically using the 4th-order Runge-Kutta 
method. 

In contrast to the constant viscosity model, for this model the viscosity varies across the channel 
width, i.e., the mean flow is characterized by a stratified viscosity. It is straightforward to verify 
that the viscosity at the lower wall increases with increasing Mach number, and hence the degree 
of viscosity stratification increases with increasing M. 



III. LINEAR STABILITY ANALYSIS 

For the linear stability analysis, the mean flow, qo = {Uq, 0, 0, po, Tq)^, is perturbed with small 
amplitude perturbations q = qo + q, and the governing equations ([U) to Q are linearized around 
the mean flow. Seeking normal mode solutions of the resulting linearized partial differential equa- 
tions, 

q(a;,2/, = q'(?/) exp [i{ax + f3z - iut)], (11) 
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we obtain a differential eigenvalue system: 

£q' = u;Iq', (12) 

where C is the linear stability operator, q' = {u', t>', w', p', T'}^ is the eigenfunction and I the 
identity matrix. Here a and /? are the stream-wise and span-wise wave-numbers, respectively, and 
= cUr + icuj is the complex frequency; the phase speed of perturbation is given by = LOyj a and 
the growth/decay rate by cjj. 

The boundary conditions on perturbation variables are taken to be: 

u'{0) = = 
v'(0) = = v'(l) 

(13) 

w'{0) = = w'{l) 
T'{1) = = f^(0). 



The Chebyshev spectral method lllSh is used to discretize the differential eigenvalue problem (fT2l - 
[T3l) at (N + 1) Gauss-Lobotto collocation points, where is the degree of the Chebyshev poly- 
nomial. This yields an algebraic eigenvalue system, AX = ujBX, which is then solved using the 
QR-algorithm of the Matlab software. 



A. Spectra and acoustic modes 

Figure 1(a) shows the distribution of eigenvalues, c = uj / a = Cr + ici, in the complex plane, 
and the zoom of Fig. 1(a), portraying the well-known 'Y' -branch of the viscous modes, is shown 
in Fig. 1(b). The parameter values are set to Re = 10^, M = 5, a = 0.1 and P = 0, with 

= 150. The classification of inviscid eigenvalues (i.e., acoustic modes) into odd- and even- 
families in Fig. 1(a) is based on their phase speeds [il9il : the odd-modes (I, III, ...) have phase 
speeds greater than unity in the limit of a 0, and the even-modes (II, IV, ...) have phase 
speeds less than zero as a; ^ 0. (Recall that the non-dimensional velocity of the top and bottom 
walls are 1 and 0, respectively.) With increasing a, however, the phase speeds of even/odd modes 
increases/decreases (not shown), and these modes become unstable once they enter the viscous 
range of the spectra (i.e., for < c,. < 1) for a range of supersonic Mach numbers and above 
some critical value of Reynolds number (see below). More specifically, the phase speed of mode I 
decreases below unity and that of mode II increases above zero, when they degenerate into unstable 
modes with increasing a. This overall scenario of modal- structure holds for both mean flows; 
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FIG. 1: Distribution of eigenvalues {c = uj/a = Cr + icj) in the complex plane for uniform shear flow with 
Re = 10^, M = 2, a = 0.1 and /? = 0. Panel (b) is the zoom of the viscous modes in panel (a). According 
to the phase-speed based classification of inviscid modes, the mode-Hi stays on the right of the mode-I, the 
mode-IV is on the left of the mode-II, and so on (see text for details). 

however, there are important differences with regard to the unstable zones in different control 
parameter space, the dominant instability and the critical Reynolds number as detailed below. 



B. Stability diagram and dominant instability 

Figures 2(a-c) show the contours of the growth rate of the least decaying mode, coidi — max(a;i), 
in the (M, a) -plane for the uniform shear flow with two-dimensional disturbances (/3 = 0) at three 
different Reynolds numbers. The flow is unstable inside the neutral stability contour (uidi — 0) 
and stable outside. With increasing Re, the size of the instability region increases and there is 

an additional instability loop in Fig. 2(c) for i?e = 5 x 10^. For a comparison, the analogue of 
Fig. 2(c) is displayed in Fig. 2(d) for the non-uniform shear flow. It is seen that the ranges of M 
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FIG. 2: (a-c) Stability maps for the uniform shear flow in the (M, a)-plane for two-dimensional {j3 = 0) 
perturbations at different Reynolds numbers: (a) Re = 10^; (b) = 2 x 10^; (c) i?e = 5 x 10^. Panel 
(d) is the analogue of panel (c) for the non-uniform shear flow at i?e = 5 x 10^. In each panel, the neutral 
contours {ui = 0) along with a few positive growth rate {uji > 0) contours are shown. 

and a, over which the flow is unstable, are much larger for the uniform shear flow. Moreover, 
the additional unstable loop at large a in Fig. 2(c) is missing in the stability diagram of the non- 
uniform shear flow in Fig. 2(d). Comparing the contours of positive growth-rates in Figs. 2(c) 
and 2(d), we find that the maximum growth-rate in the uniform shear flow can be larger by an 
order-of-magnitude. 

Figure 3(a) shows the variation of the most unstable mode with o; at a Mach number of M = 15, 
with other parameters as in Fig. 2(c). The solid line denotes the growth rate (uji) and the dashed 
line the phase speed (cph = Cr — cOr/ a). It is observed that the flow is stable for low a, but 
becomes unstable at a ^ 1.65, with the corresponding phase speed crossing zero which implies 
that this instability belongs to the mode-II [see Fig. 1(b)]; the flow becomes stable again for large 
enough a (> 6.2). (Below a < 1.65, the mode-I is the least-stable mode for which Cph > 1, and 
hence the phase-speed changes abruptly at a ^ 1.65.) Three peaks on the growth-rate curve in 
Fig. 3(a) correspond to three distinct instability loops in Fig. 2(c). It is observed that the phase 
speed changes smoothly across the kinks on the growth-rate curve for a > 1.65, implying that 
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FIG. 3: (a) Variations of the growth rate (cjj) and the phase speed (cph) of the most unstable mode with 
a for two-dimensional perturbations (/3 = 0) and M = 15; other parameters as in Fig. 2c. (b) Effects of 
Reynolds number, Re, on the growth rate of the most unstable mode with a for /? = and M = 15. 

there is no "mode-crossing" across these apparent kinks. Hence, all three unstable peaks belong 
to the same mode (see following paragraph), and, according to the above mode-classification, the 
origin of this instability is mode-II. The effect of Reynolds numbers on the most unstable mode 
is shown in Fig. 3(b), with parameter values as in Fig. 3(a). It is observed that increasing the 
value of Re from 5 x 10^ to 5 x 10^, an additional unstable peak appears on the growth-rate curve 
near o; = 11; however, the dominant instability (i.e., the mode having the maximum growth-rate 
for all a for given Re and M) still comes from the third peak [that corresponds to the uppermost 
instability-lobe in Fig. 2(c)], and this observation holds even at larger values of i?e = 5 x 10^. 
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FIG. 4: (a) Enlarged view of the first peak of Fig. 3b around a ~ 1.5. (b) Enlarged view of the third peak 
of Fig. 36 around a ~ 6. (c) Variation of the phase-speed curve corresponding to the growth-rates in panel 
(b). 

To find out whether the sharp peaks on the growth rate curves in Fig. 3 are bounded, we show 
the enlarged views of the first and third peaks [of Fig. 3(b)] in Figs. 4(a) and 4(b), respectively. 
It is clear that the growth rate varies smoothly across each peak, and the maximum growth rate at 
each peak is bounded; the sharpness of the first and third peaks in Fig. 3 is a consequence of large 
variation in growth-rate (albeit smoothly) over a small range of a. Figure 4(c) shows the phase- 
speed variation corresponding to the third-peak [i.e., Fig. 4(b)]; clearly, there is no discontinuity 
on the phase-speed curve too. [The phase-speed variation across the first peak in Fig. 4(a) is also 
smooth, not shown.] These results suggest that the instabiUty in Fig. 3 belongs to the same mode 
and the maximum growth-rate at each peak remains bounded. 
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From the zoom of the first peak, as shown in Fig. 4(a), we observe that the peak-height dimin- 
ishes with increasing Re- this is a viscous instability since it disappears in the inviscid limit. On 
the other hand, the height of the second, flatter, peak in Fig. 3(a) increases with increasing Re that 
eventually approaches the asymptotic results on the inviscid mode 11 instability of Duck et al. [|l9t| . 
The effect of Re on the third peak in Fig. 3(a) can be ascertained from its enlarged version in 
Fig. 4(b). This instability becomes stronger with increasing Re, implying that this is an inviscid 
instability too. It may be noted that this inviscid instability was not reported in Ref. [20] for the 
nonuniform shear flow. 
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FIG. 5: (a) An expanded view of the stability map in Fig. 2(c) for low Mach numbers, (b) Variations of the 
most unstable mode with a at M = 3. 



Figure 5(a) shows the zoom of the left hand comer of the stability map in Fig. 2(c). There are 
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two narrow loops of instability at M < 4. To find out the modal-origin of these two instability 
loops, we plot, in Fig. 5(b), the variations of the growth-rate (solid line) and the phase speed 
(dashed line) of the least-stable mode with a at M = 3. From the variation of the phase speed Cph, 
we find that the first unstable peak is due to the mode IT (phase speed near zero) and the second 
peak due to the mode I (phase speed near unity). Therefore, the upper "narrow" instability loop in 
Fig. 5(a) belongs to mode I and the lower loop to mode II. 
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FIG. 6: Variations of the maximum growth rate over all a, uf, and the corresponding phase speed, Cph, 
with M at i?e = 5 X 10^: (a) uniform shear; (b) non-uniform shear. The arrows in panel (b) are used to 
identify the portions of uf and Cph over M for both mode-I and mode-II. 



To find out the dominant instability mode over all a in Fig. 5(a), we plot the variation of the 
maximum growth rate 



maxa;j, 

a 



(14) 
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with M in Fig. 6(a) which increases monotonically with increasing M for the range of M shown 
(in fact, ujf decreases beyond a critical value of M , see Fig. 2). It is clear from the phase-speed 
variation in Fig. 6(a) that there is no "mode-crossing" and the mode II remains the dominant 
instability for all M. This conclusion is in contrast to the result of Hu and Zhong [|20|] (for non- 
uniform shear flow) who found that the mode I remains the dominant mode at small M and the 
mode II at moderate-to-large M, as it is evident from Fig. 6(b). For the non-uniform shear flow, 
the range of M over which the mode I remains the dominant mode increases marginally with 
Reynolds number (not shown for brevity). For example, at i?e = 5 x 10"^, the mode I is the 
dominant mode for M ~ (1.5 — 3) and the mode II for M > 3. 

The effect of three-dimensional perturbations on the least stable growth rate is shown in 
Fig. 7(a) for different span-wise wavenumber (3, with parameter values Re = 5x 10^ and M = 15; 
the zoom of the third-peak is displayed in Fig. 7(b). Comparing different growth-rate curves with 
the one for two-dimensional perturbations (P = 0), we find that there is a window of a, slightly 
beyond the third-peak, over which the three-dimensional perturbations are more unstable than their 
two-dimensional counterparts. Therefore, in general. Squire's theorem is not valid for the present 
flow configuration. This finding is in variance with the previous work [|2l|] that Squire's theorem 
holds irrespective of the value of a in the uniform shear flow of an "isothermal" compressible 
fluid. 



C. Critical Reynolds Number 

Figures 8(a) and 8(b) show the contours of the least stable growth rate in the {Re, a)-plane 
for two-dimensional disturbances (/5 = 0) with M = 3 and M = 5, respectively. The upper 
and lower instability loops in Fig. 8(a) correspond to mode I and mode II instability, respectively, 
whereas the instability loop in Fig. 8(b) arises solely from mode II. For M = 3, the flow becomes 
unstable to mode I at {Re, a) ^ (123900, 2.835), and to mode II at {Re, a) ^ (50060, 2.545). 
Therefore, the critical Reynolds number (Recr) at which the instability sets in first is determined 
by mode II in uniform shear flow - this observation holds at other values of M. A comparison of 
the values of Recr and acr between the uniform and non-uniform shear flows is given in Table I for 
different Mach numbers. It is clear that the critical Reynolds number for the uniform shear flow is 
significantly smaller than that for its non-uniform counterpart; for example, at M = 10, Recr of 
two mean flows differ by a factor of 5.6. Therefore, we conclude that the viscosity-stratification of 
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FIG. 7: (a) Variations of the growth rate {ui) of the most unstable mode with a for various values of the 
spanwise wavenumber (3 with i?e = 5 x 10^ and M = 15. (b) Zoom of panel (a) around the third peak. 

the base-flow would lead to a "delayed" transition in compressible Couetteflow in terms of modal 
instability. Another interesting observation in Table I is that the variation of Rtcr with M is non 
monotonia in the sense that the critical Reynolds number reaches a minimum at some intermediate 
value of Mach number. 
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FIG. 8: Stability maps for uniform shear flow in the (Re, a)-plane for two-dimensional (P = 0) perturba- 
tions: (a) M = 3; (b) M = 5. In each panel, the neutral contours (ui = 0) along with a few positive growth 
rate (u>i > 0) contours are shown. 





Uniform Shear 


Non-uniform Shear 


Mach Number 










M=3 


50060 


2.545 


164 900 


2.840 


M=5 


23 830 


2.130 


85 725 


2.570 


M=10 


45040 


1.870 


252700 


2.485 


M=15 


85150 


1.810 


655 850 


2.490 


TABLE I: Critical stability parameters for /3 = 
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The effect of Reynolds number on the mode 1 instability [upper loops in Fig. 5(a) and Fig. 8(a)] 



is to make it a neutral mode in the inviscid limit as is the case for non-uniform shear flow [11911 . 
This effect is similar to the first peak mode 11 instability in Fig. 4(a) where the viscosity plays 
a destabilizing role. Therefore, while the viscosity plays a dual role of destabilizing [at small a 
as in Fig. 4(a)] and stabilizing [at moderate-to-large a as in Fig. 4(b)] the mode 11 instability, it 



destabilizes the mode 1 instability. This conclusion also holds for the non-uniform shear flow [|20n . 

Even though we have presented all stability results on mode-1 and mode-11 instabilities, it may 
be noted out that the higher-order even (IV,...) and odd (111,...) inviscid modes can also become 
unstable but they remain sub-dominant with respect to mode-11 instability. 

D. Energy Analysis: Instability Mechanism 

The exponential instability can be understood by considering the rates of transfer of energy by 
the different terms in the momentum and thermal equations. For this we need to define a suitable 
norm of the perturbations which can represent the energy. We define the perturbation energy 
density as 

£{a,/3,t)= [ qt(y,t)A<q(y,t)di/, (15) 
Jo 

where the superscript f on any quantity refers to its conjugate value, and the weight matrix is 
diagonal and positive definite. Among various choices of the weight matrix Ai, we consider the 
following: 

M = diag{po, Po, Po, T^/polM\ po/7(7 - l)ToM''}, (16) 



that corresponds to the well-known Mack-norm \TI\ that has been used in many transient growth 



studies on compressible flows [Si, uM- A special property of this norm is that this definition of 
energy is free from any contribution due to the pressure related terms in the governing equations. 

Equation ([T5l) can be written for the least decaying mode, which has an exponential time de- 
pendence, as 

^0 

where the subscript 'Id' refers to 'least-decaying' mode. The rate of change of this energy with 
respect to time can be written as 



= 2%uju) exp[23(cu,,)t] c{,^\y)Mc[,M^y, (18) 
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which can be manipulated using equation (fT2l) to yield 

dSid 



dt 



iexp[253(a;zrf)t] / c(y^\y)MCc(Mdy + c.c. (19) 



Now, we decompose the total energy-transfer-rate into those coming and going through differ- 
ent physical routes. 



-Id 



exp[29(a;i,)t]5^4, (20) 



dt 

where the explicit forms of the £/s are given in the appendix. £q is the energy-transfer-rate due 
to the convection by mean flow, £i is the same from the mean flow to the perturbation, £2 is 
due to viscous dissipation, £3 is due to the thermal diffusion, and finally £4 is due to the viscous 
dissipation term in the thermal energy equation. 

Note that the above expressions involve the eigenfunction of the least-stable mode and its 
derivative. The numerical estimation of these quantities is a challenging one for the least-decaying 
mode at high Re and M with large a and p. The streamwise velocity and temperature pertur- 
bations exhibit boundarylayer like steep variations near the wall. These variations are extremely 
rapid at high a. Moreover, at high a there are also internal layers. An accurate estimation of the 
above quantities will require a highly resolved scheme to capture these steep variations. There- 
fore we used a multidomain spectral calculation, with appropriate matching conditions which can 



be found in [|23 



2M except that we have relaxed the matching of the derivative of the density 



perturbation, since the highest order of density is one in the continuity equation. A check on the 
accuracy of the results has been made by estimating the energy transfered by the pressure terms 
which must be vanishingly small by the definition of the Mack energy norm. 

Figure 9 shows the rates of different constituent energies routed via different physical processes 
at M = 5 for Re = Ax 10^. In this figure, £4 is not shown since it is negligibly small. Figure 9(a) 
shows results for 2D modes for a range of a. The sudden changes for 1.5 < a < 2 is due to a 
mode-crossing. The energy transfered from the mean-flow plays a dominant role for the onset of 
instability. The viscous dissipation and thermal diffusion plays the role of routing the energy out 
of perturbations; it is interesting to note that the thermal diffusion rate is dominant over the rate of 
viscous dissipation for 2D modes. Figure 9(b) shows these energy transfer rates for 3D modes for 
a range of P with a = 3. The main difference is that at high values of /3 the viscous dissipation 
dominates over thermal dissipation for 3D modes. This observation holds at other values of M 
and Re. 
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FIG. 9: Rates of transfer of different energies [£i, see Eq. (20)] for Re = 4 x 10^ at M = 5. Solid line, 
total energy-transfer rate; dashed Une, viscous dissipation; dotted line, from base-flow; dash-dotted line, 
thermal diffusion, (a) (3 = 0;(h) a = 3. 
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FIG. 10: Same as Fig. 9(a) but for M=3. 
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As shown in Fig. 8(a) there are distinct regions of instabilities in the (Re — Q;)-plane due to mode 
I and mode n. In order to study the characteristics special to each of these modes, we show the 
budget of energy-transfer-rates across a range of a spanning two different regions of instabilities in 
Fig. 10, with parameter values as in Fig. 8(a). Both mode I and mode II instability regions exhibit 
a qualitatively similar behavior in the shares of each physical processes except that the balancing 
involved is quantitatively different for each mode. For mode I instability, the energy transfer rate 
from the mean flow and the thermal-diffusion rate are much larger than those for mode II. 




FIG. 11: Energy transfer rates, £i, versus y for a = 2.75, ^ = and i?e = 4 x 10^ at M = 3. (a) Solid 
line, rate of total energy; dashed line, viscous dissipation; dash-dotted line, transfered by pressure terms; 
dotted line, rate of transfer from mean-flow, (b) solid line, rate of heat produced by friction, dash-dotted 
line, thermal diffusion rate. Insets in each panel show them at close to the lower wall. 

The distinction between mode I and mode II instability becomes clearer when we look at the 
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FIG. 12: Same as Fig. 11 but for a = 2.3. The subplots on the right show energy transfer rates near the 
upper wall. 

distribution of different perturbation energies across y. Figure 1 1 shows this distribution for a — 
2.75 [which belongs to mode I in Fig. 8(a)], and Fig. 12 shows the same for a. — 2.3 (which 
belongs to mode n). These figures show that the energy-transfer rate from mean-flow occurs close 
to the moving and isothermal top wall for mode I, whereas it occurs in the bulk of the flow domain 
for mode II. As one expects the viscous-dissipation is culminated near the walls. This is more at 
the lower wall for mode I, and at the upper wall for mode 11. Both modes exhibit a larger loss- 
by-conduction near the top wall. This is because the temperature gradient is set to zero at the 
lower wall via the adiabatic boundary condition. The heat generated due to viscous dissipation 
is more near the lower wall for mode I and higher near the upper wall for mode II. This is in 
accordance with the momentum loss due to viscous dissipation for both modes. Finally, apart 



21 



from the constituent energy-transfer rates found in appendix, Figs. 1 1 and 12 also show the local 
energy-transfer rates by pressure. The expression for this quantity, say, S^{t, y) is 

^5 it, y) = [pWv +vWp) + c.c. (21) 

Though this quantity does not contribute to the overall total energy-transfer rate (since this quantity 
vanishes upon integration across the channel width), it plays a role in the distribution of the same 
across the channel width. 



IV. TRANSIENT ENERGY GROWTH 

Let us write the linear stability equations in an evolution form: 

where q(?/, t; a, (3) is the inverse Fourier transform of q(x, t); the elements of the linear op- 
erator, £, are omitted for sake of brevity. In contrast to the modal linear stability analysis that 
deals with the long-time dynamics of any system via the normal-mode approach, the key idea of 
the non-modal analysis is to probe the short-time dynamics of the system in terms of perturba- 
tion energy in the parameter space where the flow is stable (such as in Fig. 2) according to the 
linear stability analysis, and investigate the potential of such stable flows to amplify the initial 
perturbation energy. 

Let a, (3] Re, M) be the maximum possible energy amplification at any time t, i.e., 

a, /3; Re, M) ^ G{t) = max (^3) 

q(o) £{a,fJ,0) 

where G{t) is optimized over all initial conditions which is computed using the singular value 
decomposition. For an efficient computation of G{t), only a selected portion of the spectra (see 
Fig. 1) is chosen [18], corresponding to the modes whose phase speeds are within the range —1 < 
ujr/a < 2 (i.e., comparable to the extremes of the mean flow velocity which varies between and 
1), and the decay rate is less than 0.5 (i.e., tOi > —0.5). With this choice of modes, the number of 
selected modes K (<< 5A^, where (A^ + 1) is the number of collocation points) can be reduced 
by a factor of 5 or more. The related details on numerical scheme are documented in our earlier 



paper Um . 
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FIG. 13: (a) Variation of the energy amplification factor, with time for parameter values of Re = 10^, 
M = 2 and a = 0. (b) Contours of Gmax in the (a, /3)-plane for Re = 10^ and M = 2. (c) Pattern of 
optimal perturbation velocities at i = in the (z, y)-plane for a = 0, /3 = 3, M = 2 and Re = 10^. 



A. Results on Energy Growth and Optimal Perturbations 



The variation of G{t, a, jS) with time for uniform shear flow is shown in Fig. 13(a) for different 
span- wise wavenumber (3, with i?e = 10^, M = 2 and a = 0; the solid and dash lines correspond 



23 



to /3 = 1 and 3, respectively. It is observed that the initial energy density can be amplified by 
a factor of 10^ or more over a time-scale of order t = O(IO^) for both /3; in the long-time limit 
(t — > oo), G{t) decays to zero since the flow is stable. Figure 13(b) shows the contours of the 
maximum amplification of energy over all time [that occurs ait = tmax such as in Fig. 13(a)] in 
the (a, /3)-plane, 

Gmaxla, /5; Re, M) = maxG'(t, a, /3; Re, M) (24) 

for Re = 10^ and M = 2. It is seen that larger energy amplification occurs for smaller values of 
streamwise wavenumber. For the dash line in Fig. 13(a), the optimal velocity patterns in the {y, z)- 
plane at t = is shown in Fig. 13(c). [The velocity pattern at t = tmax looks similar to that in 



m 



Fig. 13(c).] This represents a pure streamwise vortex which is typical of all shear flows [|4|,|5L 
The structural features of optimal patterns in compressible uniform shear flow look similar to those 
in incompressible shear flows. 

.5 



,x 10 




FIG. 14: Variations of (a) the optimal energy growth, Gopt, (b) optimal time, topt> with Mach number M 
for Re = 10^. The solid line correspond to the uniform shear mean flow, and the dashed line to non-uniform 
shear. 



The global maximum of Gmax over all combinations of wavenumber (a, (3), 



Gopt{Re,M) = supGmaxla,/?; Re,M), 



(25) 



is called as the optimal energy growth Gopt that occurs at (topt, Q^opt, /^opt)- The variations of 
Gopt and the corresponding optimal time topt^ with Mach number M are shown in Fig. 14(a-b). 
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The solid and dashed lines in each panel correspond to the uniform and non-uniform shear flow, 
respectively; Re = 10^ for these plots. Both Gopt and topt decrease monotonically with increasing 
M. The magnitude of Gopt is much larger for the uniform shear flow; the optimal time topt is 
also larger by a factor of two or more, implying that the energy growth can be sustained over a 
longer duration in uniform shear flow. These overall observations on transient energy growth hold 
at other sub-critical values of M and Re. Therefore, the uniform shear flow is more susceptible 
to sub-critical transitions than its non-uniform counterpart. As in the case of modal instability in 
Section III.C, we can conclude that the viscosity stratification along with non-uniform shear would 
also lead to a "delayed" subcritical transition in compressible Couette fiow in terms of nonmodal 
instability. 



B. Scalings of Gmax and t^ax 

In a recent paper [18], we have shown that the wellknown scaling law of incompressible shear 
flows 130, Gmax varies quadratically with the Reynolds number Re, and t^^x varies linearly with 
Re for streamwise-independent (a = 0) modes, does not hold for the non-uniform shear compress- 
ible Couette flow. To check the validity of this scaling law for the present uniform shear flow, we 
have plotted in Fig. 15(a) the variations of the rescaled energy growth ^jG{t)/ Re with rescaled 
time t/ Re for four different Reynolds number at M = 2 and (3 = 1.0; the corresponding plot for 
the non-uniform shear flow is displayed in Fig. 15(b). (Plots for different /? look similar and hence 
not shown.) It is clear that the quadratic scaling of Gmax with Re holds for the uniform shear case 
but does not hold for its non-uniform shear counterpart. 



For the nonuniform shear flow, we have argued UlSll that the following terms, associated with 
density and temperature fluctuations, in the y- and ^-momentum equations, 

£24 = -i {Toy + To^) /polM\ £34 = ^T^hM\ ^^^^ 
= -i(po, + Poi^)/po7M2, £35 = P/lM\ 
are responsible for the violation of the above quadratic scaling-law since setting them to zero, 
£24 = ^25 = = -^^35 = 0, the rescaled energy-growth curves for different Re collapses 
onto a single curve. Interestingly, for the uniform shear flow too, the above terms ^25, ^si 
and £35 remain non-zero, but the quadratic-scaling still holds. In this paper, we resolve this 
apparent contradiction via the following analysis of the linear operator in conjunction with the 
Mack transformation. 
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FIG. 15: Verification of tfie quadratic scaling law for G(€) for streamwise-independent (q = 0) perturba- 
tions with /3 = 1.0 and M = 2: (a) Uniform shear; (b) non-uniform shear. 



Let us rewrite the linear stability equation (22) as 

'di ~~ 

where 



i^^q - i£fV, 



/:(i+i)(j+3), for i = l,2;j = l,2, 

0, for i = 2, 3; j = 4, 5, 
\y,w\ and = {p, T}. 



a = £, with 



(27) 

(28) 
(29) 
(30) 



Note that the operator O' comes from y and z-momentum equations, with elements as in Eq. (22). 
Under the Mack transformation [22], {m, 0, ^, t] {Re u, 0, Re V', Re t}, Eq. (ITT]) transforms 
into 
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where £ is independent of Re and q = (u, (p, ■ip)'^. In terms of these barred-variables, an evolution 
equation for the total perturbation energy density (15) can be derived as 

dE - - - 

-^ = -i q^MCqdy-iRe'^ po4>^ C^ipdy + c.c, (32) 
dt Jo Jo 

where c.c. represents complex conjugate terms. This equation can be integrated with respect to i 
to yield, 

S{t) = E{t) + Re'^Spit), (33) 

where E{i) is the first term in Eq. (|32|) integrated with respect to i, and the second term, £p{i), 
represents the energy associated with operator If we divide by Re^ in Eq. (|3T1) . then Eq. (|33l) 
becomes independent of Re, and hence we expect the scaling of G{t) to hold. 

The above analysis is verified in Fig. 16 where the energy growth curves for different Reynolds 
numbers are seen to collapse on a single curve for the rescaled operator — > D' /Re^ in Eq. (|3T1) . 
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FIG. 16: Variation of ^/G/Re with time by rescaling the operator /Re^ in Eq. (31), with 

parameters as in Fig. 15(b). 

It is interesting to note in Fig. 15(b) that the scaling G{t) ~ Re^ holds at low Reynolds num- 
bers {Re < 10'') even for the non-uniform shear flow, and hence the contribution of to the 
perturbation energy must be negligible for such low Re. The latter statement can be confirmed if 
we explicitly compute the contribution of energy due to D'. Let us express the total energy density 

£{t) as 

£it) = Y: ete, ^"P[-^("^ - [' q;tA^£qLd, + c.c, (34) 



Re = 5000 
10000 




27 



where c^'s are the expansion coefficients of q : 

q(l/> t) = ^Ck exp[-iuJkt]<^M^ (35) 

k 

which can be evaluated by the singular value decomposition of the propagator of q such that 
^^(^max) = Cmax- In Eq. (l34l) . the eigenfunction q' is normalized (to make the initial total energy 
£^(0) = 1) with respect to the weight matrix A^, such that ||A^qj^|| = 1, where Ai is given by 
Ai^Ai = Ai. It is straightforward to verify from Eq. (|34|) that the contribution of the terms in 
Eq. (|26l) to the total energy is: 

X [-iv'i^Dp', + f3w[^p'k]dy + C.C. (36) 
Jo 

Figure 17(a) shows the variation of £p with time at a Reynolds number Re = 10^; the symbols, 
circle and triangle, correspond to times at which Gmax occurs for non-uniform and uniform shear 
flows, respectively. It is seen that for the case of non-uniform shear SpStt = t^aax is much larger in 
comparison with that for uniform shear. At a low Reynolds number Re = 10^, however, <£^p(tmax) 
is negligible for both uniform and non-uniform shear flows [see Fig. 17(b)], and hence the scaling 
of G{t) holds for relatively small Re [see Fig. 15(b)] in non-uniform shear flow. 

The above analysis suggests that the (streamwise-independent) linear operator C of compress- 
ible flows can be partitioned into a Reynolds number dependent operator, Cp, and a Reynolds 
number independent operator C [Eqs. (|27l) and (|3T])1. The contribution of this i?e-dependent oper- 
ator, Cp, to perturbation energy would decide whether the scaling Gmax ~ would hold or not 
for a given mean flow. For the uniform shear flow, Cp has negligible contribution to the energy 
growth and hence the quadratic scaling-law holds. 

C. Inviscid Algebraic Growth and Optimal Perturbation 

The purely inviscid nature of the algebraic growth suggests one could try to obtain the transient 
growth characteristics directly from inviscid equations. As it has been shown numerically in the 
earlier section that the algebraic growth is very pronounced for the modes that are independent of 



the streamwise coordinate (i.e., a = 0). For such an unidirectional flow, EUingsen & Palm n25\\ 



had found an analytical solution for incompressible flows. An extension of this solution for density 
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FIG. 17: Variation of £p with time for uniform and non-uniform shear flows: M = 2, /3 = 1 and a = 0. 

(a) Re = 10^ and (b) Re = lOl 



and temperature perturbations was considered for the compressible situation 12611 which resulted 
in a constraint due to the continuity equation which relates spanwise velocity with normal velocity. 
Further, this solution would also result in another constraint which relates density and temperature 
perturbations, which was not considered before, but is considered here [|24ll . The EUingsen-Palm 
solution for compressible shear flow can be written as 



^ivs 
^ivs 

Pivs 

T 



I 

^ivs 



T' 



(37) 
(38) 
(39) 
(40) 
(41) 



where u'-^^, v{^^ and 7]'^^ are the initial perturbation quantities which are to be determined via an 
optimization procedure; in the following analysis, the subscript "ivs", which refers to "inviscid 
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solution", is dropped for the sake of simplicity. The perturbation energy S{t) can be written in 
the basis of the quantities u', v' and T', after removing w' and p' using the above mentioned 
constraints, as 

m = (po\u\' + ^[po{f3' - D') - poyD]v 

^° ^fAdy. (42) 



(7 - 1)M2 

Let i/j = {u, V, T}^ and -0' = {u', v', T'}^. Then the above equation can be written as 

S{t)= [ i^'^A^MAiP'dy, (43) 
Jo 

where M = diag{po, [poiP'^ — D^) — PoyD]//3'^, pl/{j — 1)M^}, and A is a 3 x 3 matrix which can 
be defined by casting Eqs. (l37l) . (l38l) and (|4T]) in the form, -0 = A?/''. Now G{t) = max^' is 
given by 

=max({Aj), (44) 
where A^'s are the eigenvalues of the differential equation 

A^MA^;' = \Mii)' (45) 



with the boundary conditions t>'(0) = f '(1) = 0. In contrast to Hanifi & Henningson's [|26l] four- 
variable model, this equation (|45] ) has only three dependent variables and hence called a "reduced" 
model. The constraint of vanishing pressure fluctuation is essential to obtain this reduced model; 



the related spatial problem has been solved elsewhere 112411 . 

Eq. (l45l) has been solved using the spectral method. Figure 18 shows the inviscid algebraic 
growth curve G{t) at M = 5 and /3 = 1. The viscous transient growth curves are also shown for 
three different Reynolds numbers. It is seen that for the entire growth duration the viscous and 
inviscid growths agree quantitatively, demonstrating the inviscid nature of the algebraic growth. 



In terms of energy-transfer-rate, only the following term (see appendix): 

m = - 



^ TT ~t~ , ^OPO?/ ~t~ 

PoUoyU^V H ITj^P^ 







ro7(7 - 1)M2 



PolM^' 

dy + c.c. (46) 



survives in the inviscid limit. It is clear that the energy transfer from the mean flow occurs via the 
Reynolds stress (u'^v) and the coupling of the normal perturbation velocity with density (p'^v) and 
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FIG. 18: Viscous and inviscid energy growth curves for /3 = 1 and M = 5. Solid line, inviscid solution, 
G{t); dashed line, full viscous solution, G{t). 

temperature (T^v). The last two contributions {p^v and T^v) are unique to compressible flows. 
Further, Eqs. (|37] ). (|40l ) and (|4TI) also suggest that this inviscid growth is due to the transfer of 
energy from mean flow to u, p and T via the fluctuation in the normal velocity, v. The continuity 
is satisfied by a mere readjustment of w which need not grow due to this algebraic growth. The 
growth of u eventually would give rise to streaks. 

Figure 19 shows the optimal patterns of the perturbation velocity-field at t = 100, obtained 
from our reduced inviscid model, Eq. (45). Figure 19(a) shows the counter-rotating stream-wise 
vortices in the (y — z)-plane, and Fig. 19(b) shows the contours of streamwise velocity fluctuation u 
in the same plane which exhibits the well-known streaks. The structural features of these vortices 
and streaks are strikingly similar to those obtained from the solution of full viscous equations. 
Therefore, the compressible inviscid EUingsen-Palm solution, along with the constraint of null 
pressure fluctuations, captures all essential features of the algebraic growth of the full viscous 
equations. 

V. SUMMARY AND CONCLUSION 

The compressible plane Couette flow is linearly unstable to acoustic disturbances for a range 
of supersonic Mach numbers. We found that the effects of viscosity-stratification and nonuniform 
shear rate are to stabilize the unstable modes at large stream-wise wavenumber (a) and Mach 
numbers (M). For a given Mach number, the critical Reynolds number {Re) is found to be signif- 
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FIG. 19: Inviscid optimal patterns of velocity fluctuations given by Eq. (l45l) for M = 5 and [3 = 1. (a) 
Velocities my — z plane; (b) contours of u. 

icantly smaller (by a factor of 3 or more) in the uniform shear flow than in its nonuniform shear 
counterpart; for a given Re, the maximum growth rate (over all a, cf. Fig. 2) could be larger 
by an order-of-magnitude in the former. This strong stabilization effect appears to be tied to the 
strong viscosity stratification in non-uniform shear flow, and, therefore, the viscosity stratifica- 
tion would lead to a delayed transition in compressible Couette fiow. Three-dimensional modes 
could be more unstable than their two-dimensional counterparts for some values of a, and hence 
Squire's theorem is, in general, not valid for the "nonisothermal" compressible Couette flow. It 
is shown that the mode n remains the dominant instability (i.e., the mode having the maximum 
growth-rate over all a, Eq. 14) for all Mach numbers in the uniform shear flow. In contrast, for 
the nonuniform shear flow, the mode I is the dominant instability for low Mach numbers and the 
mode II for moderate-to-large Mach numbers. For both mean flows, the viscosity plays the dual 
role of destabilizing (at small a) and stabilizing (at moderate-to-large a) the mode II instability, 
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but it destabilizes the mode-I instability. The higher-order odd (III,...) and even (IV,...) inviscid 
modes could also become unstable, but they remain subdominant with respect to mode I and mode 
II instabilities. 

An analysis based on the perturbation kinetic energies transfered by different terms of the gov- 
erning equation has been carried out to understand the origin of modal instabilities. The instability 
is primarily caused by an excess transfer of energy from mean-flow to perturbations for a band 
of stream- wise wavenumbers. It is found that the energy-transfer rate from the mean-flow occurs 
close to the moving and isothermal top- wall for mode I, whereas it occurs in the bulk of the flow 
domain for mode II. For 2D modes, the thermal-diffusion process tends to stabilize the fluctua- 
tions at a higher rate than the viscous dissipation; for 3D modes, however, the viscous dissipation 
dominates over thermal diffusion at high spanwise wavenumber. 

For the transient growth analysis, it is shown that the maximum temporal growth of perturbation 
energy, Gmax^ and the corresponding time scale to attain this maximum, tmax^ are much larger (and 
can differ by a factor of 5 or more) for the uniform shear flow in comparison with the nonuniform 
shear flow. (In other words, the viscosity stratification has a strong stabilizing effect on transient 
energy growth.) Therefore, the uniform shear flow is more susceptible to subcritical transitions 
than its nonuniform shear counterpart. For both mean flows, the optimal energy growth, Gopt^ (i-e-> 
the global maximum of Gmax in the (a, /3)-plane for given Re and M) decreases with increasing 
M; pure streamwise vortices («opt = 0) are the optimal velocity patterns at large M but the 
modulated streamwise vortices (ttopt 7^ 0) are optimal patterns for low-to-moderate values of 
M. The physical mechanism of transient energy growth is tied to the transfer of energy from the 
mean flow to perturbations via the Reynolds stress and the coupling of density and temperature 
perturbations with the normal velocity. 

For the streamwise independent perturbations (a = 0), we have found that the transient energy 
growth follows the well-known scaling law, Gmax ~ R^^ and tmax ~ Re-, of incompressible shear 
flow [3]. This is in stark contrast to the result on the nonuniform shear flow for which the above 
scaling law does not hold uM- An analysis of the linear stability operator, £, shows that C can 
be partitioned into a i?e-dependent operator, Cp, and a i?e-independent operator, C, (Eqs. 27 and 
31) via the Mack transformation. The (in)validity of the above scaling laws for the (non-)uniform 
shear flow is shown to be tied to the (non-)negligible contribution (to perturbation energy) of Cp. 
Lastly, a 'reduced' inviscid model (Eq. 45), based on the inviscid EUingsen-Palm-type solution, 
has been derived which captures all salient features of transient energy growth of full viscous 
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equations. 



APPENDIX A: EVOLUTION EQUATION OF PERTURBATION ENERGY 



It can be verified that the perturbation energy £{a,/3, t) satisfies the following time-evolution 



equation [|18|1 : 

d£_ 
'dt 



PoToy 



ro7(7 - 1)M2 



dy + c.c. 



£2 



Re 



\a^{pQ + \q)u^u + po{a^ + 0^)u^u 



—u\pQyD + pqD'^)u — iau^poy 

+ {po + ^o)D)v + a(3{pQ + \q)v)w 

-{UoyyPT + UoyToyPTT^T - UoyPru'^DT 

-iav^Xoy + (/io + Xo)D)u 

+Po{a'^ + P^)v'^v - v'^iiXoy + Poy)D 

+ (Ao + Po)D^ + PoyD + poD^)v 

— i/5(Ao + po)v'^ Dw — iaUoyPrv'^T 

—ipXoyv'^w — i(3pQyw'^v 

-\-a(3{pQ + Xq)w'^u — i/5(Ao + po)w^ Dv 

+ (/io(a^ + /?^) + /?^(Ao + Pq))w^w 

—Pqw'^D'^w — pQyw'^ Dvo\ dy + c.c. 

+ 2ToyPTD - {a^ + P^)po + poD^] fdy + c.c. 



(Al) 



-i / q'^A^^q dy + c.c. = £0 + £1 + £2 + £3 + ^4, 

The constituent energy transfer rates, £^-£4, have following forms (with D = d/dy): 

£0 = -ia / Uoq^Mqdy + c.c. (A2) 
Jo 



(A3) 



(A4) 



(A5) 
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2iioUoyf^Du 
+ 2iaijLoUoyT^v + U^/jLtT^T 



dy + c.c. 



(A6) 



Here, Si is the energy transfer rate from the mean flow, S2 the viscous dissipation rate, S3 the 
thermal diffusion rate and S4 the shear- work rate, respectively; note that the convective transfer of 
perturbation energy by the mean flow, Sq, is zero. 
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